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ABSTRACT 

Regional droughts are common in North America, but pan-continental droughts extending across multiple 
regions, including the 2012 event, are rare relative to single-region events. Here, the tree-ring-derived North 
American Drought Atlas is used to investigate drought variability in four regions over the last millennium, 
focusing on pan-continental droughts. During the Medieval Climate Anomaly (MCA), the central plains 
(CP), Southwest (SW), and Southeast (SE) regions experienced drier conditions and increased occurrence of 
droughts and the Northwest (NW) experienced several extended pluvials. Enhanced MCA aridity in the SW 
and CP manifested as multidecadal megadroughts. Notably, megadroughts in these regions differed in their 
timing and persistence, suggesting that they represent regional events influenced by local dynamics rather 
than a unified, continental-scale phenomena. There is no trend in pan-continental drought occurrence, de- 
fined as synchronous droughts in three or more regions. SW, CP, andSE (SW+CP+SE) droughts are the most 
common, occurring in 12% of all years and peaking in prevalence during the twelfth and thirteenth centuries; 
patterns involving three other regions occur in about 8% of years. Positive values of the Southern Oscillation 
index (La Nina conditions) are linked to SW, CP, and SE (SW+CP+SE) droughts and SW, CP, and NW 
(SW+CP+NW) droughts, whereas CP, NW, and SE (CP+NW+SE) droughts are associated with positive 
values of the Pacific decadal oscillation and Atlantic multidecadal oscillation. While relatively rare, pan- 
continental droughts are present in the paleo record and are linked to defined modes of climate variability, 
implying the potential for seasonal predictability. Assuming stable drought teleconnections, these events will 
remain an important feature of future North American hydroclimate, possibly increasing in their severity in 
step with other expected hydroclimate responses to increased greenhouse gas forcing. 


1. Introduction 

During the summer of 2012, over half of the conti- 
nental United States experienced moderate-to-severe 
drought conditions, with the core of the drought area 
extending across Mexico and the desert Southwest, the 
central plains, and the southeastern United States (Fig. 1). 
At its peak in July, the drought covered 61.8% of the 
contiguous United States (NCDC 2013a), exceeding the 
99th percentile of drought-covered area observed during 
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the full period of observation (from 1895 to the present). 
Notably, the area of the United States classified by the 
National Climatic Data Center (NCDC) as “moderately 
to extremely dry” exceeded the July 2012 area in only 
eight other months during the instrumental record (see 
http://wwwl.ncdc.noaa.gov/pub/data/cmb/sotc/drought/ 
2012/13/uspctarea-wetdry-mod.txt): April-October 1934 
and December 1939, all during the Dust Bowl drought of 
the 1930s. 

Drought is a common and recurrent feature of North 
American climate (e.g.. Cook et al. 2007; McCabe et al. 
2004; Nigam et al. 2011; Schubert et al. 2004a, b; Seager 
et al. 2005), but pan-continental droughts like the drought 
of 2012 are rare relative to single-region events. Despite 
their rarity, however, these events present significant 
management challenges because they range across broad 
regions and communities with different water resource 
constraints (e.g., irrigation from rivers versus groundwater). 
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Fig. 1. Model-estimated JJA standardized soil moisture anom- 
alies for 2012, from the CPC (Fan and van den Dool 2004). Base 
period for the standardization was 1948-2012. Analyzed drought 
regions in this study are outlined by the black dashed lines. The 
regions are defined as the southwest (SW; 32‘M0°N, 125°-105°W), 
central plains (CP; 34°M6°N, 102°-92°W), northwest (NW; 42°- 
50°N, 125°-110°W), and southeast (SE; 30°-39°N, 92°-75°W). 


ecosystems (e.g., forests and grasslands), and crops. The 
spatially extensive nature of pan-continental droughts 
also makes understanding their dynamics difficult. By 
definition, they encompass geographically broad regions 
of North America with often widely different climates, 
seasonalities, and climate mode teleconnections. Droughts 
in the Southwest, for example, are typically caused by 
winter and early spring precipitation deficits forced by 
cold sea surface temperature (SST) anomalies in the 
eastern tropical Pacific (e.g., Herweijer el al. 2006; 
Hoerling et al. 2009; Schubert et al. 2009; Seager et al. 
2005). The Northwest has similarly strong Pacific SST 
teleconnections, but of opposite sign, resulting in hy- 
droclimate anomalies that are often out of phase with 
the Southwest (McCabe et al. 2004, 2008; Steinman et al. 
2012). The central plains region, while still influenced to 
some degree by the Pacific (Schubert et al. 2004a, b), is 
on the northern edge of the Pacific SST teleconnection, 
and droughts in this area are driven primarily by spring 
and summer precipitation reductions often linked to 
warm Atlantic SSTs (e.g., Kushnir et al. 2010; Nigam 
et al. 2011). Much like the central plains, droughts in the 
southeastern United States also respond to Atlantic and 
Pacific SST variations (Mo and Schemm 2008; Pederson 
et al. 2012; Seager et al. 2009). 

Proxy reconstructions of drought allow studies of 
drought variability over much longer time scales than 
are available from the limited observational record. 
Improved sampling and characterization of multidecadal 
hydroclimate variability is therefore possible, including 
the occurrence of rare and extreme events such as 


pan-continental droughts. Additionally, paleo perspec- 
tives provide glimpses into different hydrological-mean 
states such as the often discussed megadroughts of the 
medieval era, periods of persistent (decadal and longer) 
drought that have few comparable analogs during the 
last 500 years (Cook et al. 2004, 2007, 2010; Gray et al. 
2011; Herweijer et al. 2007; Meko et al. 2007; Routson 
et al. 2011; Stahle et al. 2007; Stine 1994; Woodhouse 
and Overpeck 1998; Woodhouse et al. 2010). Herein, we 
use a tree-ring-derived, spatially explicit, proxy recon- 
struction of North American hydroclimate variability to 
investigate the occurrence and forcing of pan-continental 
droughts over the last millennium. We focus our analy- 
ses on three principal research questions: 1) How often 
have pan-continental droughts occurred over the last 
1000 years? 2) How variable is the occurrence of pan- 
continental droughts, and are certain spatial patterns 
more prevalent than others? 3) How strongly can the 
occurrence of these pan-continental droughts be linked 
to known teleconnections and modes of SST forcing? 

2. Methods and data 

a. The North American Drought Atlas 

We use an updated version of the North American 
Drought Atlas (NAD A; Cook et al. 1999, 2004), a tree- 
ring reconstruction of the summer season [June- August 
(JJA)] Palmer drought severity index (PDSI; Palmer 
1965). The PDSI is a locally normalized index of drought, 
reflecting the balance between moisture supply (precipi- 
tation) and demand (evapotranspiration as a function of 
temperature). Positive values of PDSI indicate wetter- 
than-normal conditions and negative values indicate 
drier-than-normal conditions. Because PDSI has a mem- 
ory time scale of 12-18 months (Guttman 1998; Vicente- 
Serrano et al. 2010), values during the summer season will 
reflect temperature and precipitation anomalies from 
throughout the year (e.g., St. George et al. 2010). The 
employed version of the NADA has improved areal 
coverage and spatial resolution (0.5°), allowing for more 
precise regional characterizations of drought variability. 
To achieve this, the tree-ring network used for recon- 
struction was increased by over 1000 tree-ring chronol- 
ogies, many of which were from several hundred to over 
a thousand years long. In every case, the tree-ring 
chronologies were detrended using methods that ex- 
plicitly preserved centennial-or-longer time scale vari- 
ability due to climate. The success in doing this is evident 
in Cook et al. (2010), where two centennial-duration 
megadroughts in California, discovered by Stine (1994) 
using totally independent methods, were captured with 
great fidelity. The point-by-point method of PDSI 
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Fig. 2. Spearman’s rank correlation p between the regional time series and PDSI at every grid point in the 
0.5°-resolution NADA for 1000-2005 CE in the (a) SW, (b) CP, (c) NW, and (d) SE. 


reconstruction, which uses a search radius of 450 km, 
also guarantees that there is little overlap in the tree- 
ring data used to reconstruct gridded drought in the 
four selected regions analyzed here (Fig. 2), thus 
making their comparisons with forcings at both high 
and low frequencies largely independent. While data 
from the NADA are available for the past 2000 years, 
we restrict our analyses to 1000-2005 Common Era 
(CE) to ensure consistent spatial and temporal cov- 
erage over our regions of interest. For the purposes of 
some temporal comparisons, we split the data into two 
eras: the Medieval Climate Anomaly (MCA; 1000- 
1500 CE) and the Little Ice Age/modern era (LIA/ 
MOD; 1501-2005 CE). The years bounding these two 
periods are loosely defined in the literature, and global 
analyses generally define the terminus of the MCA as 
some time during the thirteenth or fourteenth centuries 
(e.g., Graham et al. 2011; Mann et al. 2009). In North 
America, however, the megadroughts were the de- 
fining characteristic of medieval climate, and many 
of these events extended well into the fifteenth 
century. We therefore define the MCA in our anal- 
ysis as 1000-1500 CE, and the LIA/MOD period as 
1501-2005 CE. 


b. Climate indices 

To diagnose possible causes of pan-continental droughts 
in the NADA, we use three climate indices represen- 
tative of the major SST modes linked to North American 
droughts. The Southern Oscillation index (SOI; Ropelewski 
and Jones 1987; Trenberth and Hoar 1996) is defined as 
the normalized difference in sea level pressure between 
the island of Tahiti and the city of Darwin in Australia, 
reflecting variations in the state of the El Nino-Southern 
Oscillation (ENSO) system. Positive values indicate La 
Nina conditions, typically associated with drought in the 
Southwest, southern plains, and the Southeast. We use 
the version provided by the Climate Prediction Center 
(CPC) of the National Weather Service, available monthly 
from 1882 to the present. As a test for how sensitive our 
results were to our choice of ENSO indicator (the SOI), 
we repeated all teleconnection analyses with the Nino-3. 4 
index, an SST-based indicator of ENSO variability. 
Analyses with the Nino-3.4 index yielded nearly iden- 
tical results to our analyses based on the SOL The Pacific 
decadal oscillation (PDO; Mantua et al. 1997; Mantua 
and Hare 2002) is defined as the leading principal com- 
ponent of monthly SST variability poleward of 20°N in 
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the North Pacific Ocean. Positive phases of the PDO are 
generally linked to drought in the Northwest and wet 
conditions in the Southwest (McCabe et al. 2004, 2008). 
The degree to which the PDO is separable and distinct 
from ENSO is poorly understood, and there is some 
evidence that it may be simply a low-frequency, extra- 
tropical expression of ENSO (McCabe et al. 2008; 
Newman et al. 2003). While the physical mechanisms 
involved are therefore not fully elucidated, the impacts 
of the PDO on drought variability in North America 
are nevertheless well established (e.g., McCabe and 
Dettinger 2002; McCabe et al. 2004, 2008). We use the 
version provided by the Joint Institute for the Study of 
the Atmosphere and Ocean (University of Washington), 
available monthly from 1900 to the present. The At- 
lantic multidecadal oscillation (AMO) index is calcu- 
lated as Atlantic SSTs averaged over 0°-60°N, 0°-80°W 
minus the global rise in SSTs averaged over 60°S-60°N 
so that it represents internal variability independent 
of any long-term greenhouse gas-forced trend (after 
Trenberth and Shea 2006). Positive values of the AMO, 
indicative of warm Atlantic SSTs, are typically asso- 
ciated with dry conditions in the central and southern 
plains and the Southeast (Kushnir et al. 2010; McCabe 
et al. 2004, 2008; Nigam et al. 2011). This version of 
the AMO index is calculated using the National 
Oceanic and Atmospheric Administration (NOAA) 
extended reconstructed SST dataset (version 3b), avail- 
able monthly from 1854 to the present (Smith and 
Reynolds 2003). 

c. Analysis 

Within the conterminous United States, hydro- 
climatic variability is a recurrent, robust feature of the 
regional climate in four primary regions: the Southwest 
(Hoerling et al. 2009; Woodhouse et al. 2010), the 
Northwest (Nelson et al. 2011; Steinman et al. 2012), the 
central plains (Hoerling et al. 2009; Schubert et al. 
2004a, b), and the Southeast (Pederson et al. 2012; 
Seager et al. 2009). For the subsequent analyses of pan- 
continental drought, we spatially average the NADA 
over these four regions to derive regional time series of 
PDSI (1000-2005 CE). Regional boundaries were cho- 
sen to maximize statistical independence between the 
PDSI time series, while still encompassing the major 
areas of drought variability: SW: 32°-40°N, 125°-105°W; 
CP: 34°-46°N, 102°-92°W; NW: 42°-50°N, 125°-110°W; 
and SE: 30°-39°N, 92°-75°W. We do not expect these 
four regions to be completely independent because of 
various processes that govern the spatial expression of 
drought, including topography, regional circulation, and 
various teleconnection modes. The correlation maps (p) 
between these regional-average time series and PDSI at 


each grid point in the NADA (Fig. 2) nevertheless 
strongly indicate that positive correlations are concen- 
trated within the bounding boxes of the defined regions 
(indicated by the black dashed boxes) and decay rapidly 
outside of them. These time series therefore can be in- 
terpreted as reasonably distinct and separable regions of 
drought variability, appropriate for our analysis. Droughts 
are characterized to have occurred in the regional-mean 
time series when PDSI falls to a value of —0.5 or lower in 
any individual year (a value of —0.5 being considered an 
incipient drought; see definition at http://www.esrl.noaa. 
go v/psd/data/usclimdi vs/data/). We used a modest thresh- 
old of —0.5 rather than a more extreme value (e.g., —1.0) 
to ensure adequate sampling of events over the last 
millennium and during the overlapping period with 
the teleconnection indices. A PDSI threshold of —0.5, 
for example, yields 124 SW, CP, and SE (hereafter, 
SW+CP+SE) pan-continental drought events over the 
last 1000 years, while a threshold of —1.0 yields only 47. 
During the instrumental period (1850-2005), a shift 
from —0.5 to —1.0 changes the number of events from 23 
to 5, drastically reducing our ability to sample the cli- 
mate indices. With these considerations in mind, and 
given that NOAA refers to PDSI values of —0.5 as an 
incipient drought, we feel that a PDSI value of —0.5 is 
a statistically reasonable and physically defensible thresh- 
old. We define a pan-continental drought as any in- 
stance when any three [SW+CP+SE; SW, CP, and NW 
(hereafter, SW+CP+NW); SW, NW, and SE (hereaf- 
ter, SW+NW+SE); or CP, NW, and SE (hereafter, 
CP+NW+SE)] or all four [SW, CP, NW, and SE 
(hereafter, SW+CP+NW+SE)] of our defined regions 
simultaneously have PDSI values of —0.5 or lower in the 
same year. By this definition, the four-region droughts 
will overlap with and also be counted as three-region 
droughts. While it would be possible to exclude the four- 
region events from the three-region drought counts, we 
have chosen to allow this overlap to maximize the sam- 
pling of these relatively rare events. This sampling is es- 
pecially important for the teleconnection analysis, which 
is limited to the observational period. 

For both the single-region and pan-continental droughts, 
we analyze temporal changes in the drought occurrence 
and average aridity in the NADA, and use the obser- 
vational climate indices (SOI, PDO, and AMO) to es- 
tablish possible forcing mechanisms. We correlate (p) 
the dynamic indices against each NADA grid point 
during their periods of overlap, focusing on the main 
seasons of influence associated with each climate mode: 
December-February (DJF) and March-May (MAM) 
for the SOI and PDO; and MAM and JJA for the AMO. 
Because the PDO and AMO represent low-frequency 
modes of variability (decadal and longer), the PDSI and 
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climate indices are smoothed with a 10-yr locally 
weighted scatterplot smoothing (LOWESS) filter before 
conducting these correlations. 

It is possible that pan-continental droughts do not 
arise primarily from large-scale forcing that organizes 
anomalies over large spatial scales (e.g., the AMO), but 
rather through the random superposition of drought 
events occurring independently in the four regions. 
To test this, we conduct a 5000-member ensemble re- 
sampling of each of the regional PDSI time series. In 
each ensemble member, we randomize the phase in- 
formation of the regional PDSI time series [using the 
method of Ebisuzaki (1997)], generating new series in 
each region with identical lengths and power spectra, but 
different relative phasing of their spectral components. 
Any dependencies across regions are therefore removed 
and each new time series can be interpreted as an in- 
dependent random draw from an underlying parent 
population defined by the autocorrelation structure of the 
original regional time series. We calculate the percentage 
occurrence of pan-continental droughts from these syn- 
thetic series, and then calculate the 95th percentile of 
occurrences across the entire ensemble. If the observed 
percent occurrence of pan-continental droughts in the 
NADA exceeds this 95th percentile threshold, we reject 
the null hypothesis that pan-continental droughts arise 
randomly from independent variability in each region. 

To determine the significance of the drought telecon- 
nections (SOI, PDO, and AMO), composite averages 
are computed for those climate index values associated 
with the single-region and pan-continental droughts. A 
5000-member ensemble resampling of the climate in- 
dices is then performed to generate 90th and 95th per- 
centile confidence limits. For example, during the full 
overlapping time period between the SOI data and the 
NADA (1882-2005), there are 41 years that qualify as 
SW droughts (Table 1). An average of the SOI values for 
these 41 years gives a composite SOI anomaly associ- 
ated with these events. We then draw 41 random years 
from the SOI time series and average them, repeating 
this process 5000 times. If the original composited SOI 
anomaly exceeds the 90th or 95th percentile thresholds 
of the ensemble resampling, the drought-teleconnection 
association is characterized as marginally or fully sig- 
nificant, respectively. 

3. Results 

a. Single-region drought comparisons 

The CP and SW are the prime epicenters for North 
American megadroughts, and long-term droughts dur- 
ing the MCA in these regions can be clearly seen in the 


Table 1. The number of single-region and pan-continental 
drought years (PDSI s -0.5) in the NADA available from the 
overlapping periods with the climate indices from the observa- 
tional era. The climate indices are current up through 2012, but the 
NADA ends in 2005. Total overlap is 124 years for the SOI, 106 
years for the PDO, and 152 years for the AMO. 



SOI 

PDO 

AMO 


Single-region droughts 


SW 

41 

36 

57 

CP 

34 

31 

53 

NW 

41 

36 

52 

SE 

43 

37 

59 

Pan-continental droughts 




SW+CP+SE 

16 

13 

23 

SW+CP+NW 

8 

7 

12 

SW+NW+SE 

8 

7 

11 

CP+NW+SE 

12 

11 

17 

SW+CP+NW+SE 

7 

6 

9 


regional time series (Figs. 3a, b). Megadrought activity 
in the SW and CP peaks during the exceptionally dry 
twelfth, thirteenth, and fifteenth centuries. Notably, meg- 
adroughts across these two regions are not completely 
synchronous in time, and droughts in the CP have much 
longer persistence time scales than in the SW. This 
suggests that the megadroughts do not represent a single, 
unified phenomenon dominating western hydroclimate 
during the MCA. Rather, they behave as regionally 
distinct phenomena governed primarily by particular 
teleconnections and local land surface feedbacks, pos- 
sibly with some longer time scale forcing making them 
more likely to occur during the MCA than during the 
more recent LIA/MOD period. Both regions transition 
to wetter mean conditions in more recent centuries. 

Megadroughts are generally absent in the PDSI anom- 
alies for the NW (Fig. 3c). Multidecadal pluvials are 
instead apparent during ca. 1090-1120, 1260-1350, and 
1440-1475 CE. There is additionally little change in NW 
drought variability across the transition from the MCA 
to the LIA/MOD. In the SE (Fig. 3d), interannual 
drought variability is muted compared to the SW and 
CP, and conditions are dry (in a regionally normalized 
sense) up to the middle of the nineteenth century, at 
which point a wetting trend to the present day begins. 
Independent corroboration of this long-term trend to- 
ward more pluvial conditions is difficult because few, if 
any, reconstructions of drought from this region are 
available that use completely independent proxy data 
from the NADA. However, this trend is apparent in 
more recent regional reconstructions that incorporate 
additional proxies independent of the NADA (Pederson 
et al. 2012, 2013), and the trend is consistent with in- 
dependent studies of forest dynamics over the last sev- 
eral hundred years (McEwan et al. 2011). 
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Fig. 3. Area-averaged PDSI from the NADA for the (a) SW, (b) CP, (c) NW, and (d) SE 
regions, as shown in Figs. 1 and 2. Green and brown bars are the original data, and dark black 
lines are a smoothed version of the time series using a 10-yr LOWESS spline. 


Three of the regions (SW, CP, and SE) show long- 
term, multicentennial trends in drought frequency or 
occurrence (Fig. 4a), calculated from the eleventh to 
twentieth centuries. The number of droughts per century 
has declined significantly ( p < 0.01, based on best-fit 


linear least squares regression) in the CP and SE, on the 
order of —2.2 and — 1.5 century -1 , respectively. There 
are also declines in the SW (— l.Ocentury -1 ), although 
this trend is not significant ( p = 0.11). In step with the 
declines in drought frequency, mean PDSI has been 
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No. Droughts (PDSI s-0.5) Mean PDSI 



Fig. 4. Drought statistics derived from each of the regional time series: (a) number of drought years (PDSI £ -0.5) 
per century and (b) mean PDSI calculated across all years of each century. 


increasing for the same three regions, indicating a 
trend toward wetter average conditions (Fig. 4b). The 
CP shows the strongest trend (+0.10 century , p < 
0.01), while trends in the SW (+0.06 century -1 ) and SE 
(+0.04 century -1 ) are weaker and only marginally sig- 
nificant ( p < 0.10). Trends in the NW are small and in- 
significant for both drought frequency and mean PDSI. 

The extended periods of drought in the SW, CP, and 
SE during the MCA (Fig. 2) may indicate enhanced 
drought persistence, possibly due to sustained Atlantic 
or Pacific SST forcing (Burgman et al. 2010; Conroy 
et al. 2009; Feng et al. 2008; Oglesby et al. 2012; Seager 
et al. 2007a, 2008) or land surface feedbacks (Cook et al. 
2013). We calculated the autocorrelation function (ACF) 
for each region separately for the MCA and LIA/MOD 
periods, after first adjusting each PDSI time series to 
a mean of zero (Fig. 5). Of all four regions during the 
MCA, the CP has the largest lag-1 persistence (Fig. 5b), 
twice as high as the SW, the other region of mega- 
drought activity. Persistence during the MCA in the CP 
remains consistently high and generally significant out to 
6 years, reflecting the extended, multidecadal periods of 
drought that defined climate in this region during the 
MCA (Fig. 2b). The SW has no significant autocorrela- 
tion between years 2 and 6, but does have significant 
autocorrelation at lags of 7 and 8 years, possibly due to 
the strong decadal and cyclic drought variability in this 
region. Notably, the SE does not have any significant 
persistence during the MCA beyond 1 year. Drier con- 
ditions in the SE during the MCA therefore represent 
the beginning of a long-term wetting trend, rather than 
a fundamental change in the underlying variability. 
During the last 500 years, significant persistence beyond 
a lag of 1 year disappears for all regions, although 
droughts in the CP continue to have the highest lag-1 
autocorrelation of any region. 


b. Incidence of pan-continental droughts 

The percent of total drought years for each of the five 
pan-continental drought configurations during 1000-2005 
CE are indicated in the legend of Fig. 6. The single 
most common pan-continental drought pattern over 
the last thousand years is SW+CP+SE, a pattern quite 
similar to the 2012 drought (see Fig. 1). The second most 
common pattern is CP+NW+SE; followed by SW, CP, 
and NW (hereafter, SW+CP+NW); and SW+NW+SE. 
The incidence of drought in all four regions simultaneously 
is only 5.2% of years. There is substantial centennial-scale 
variability in the occurrence of these droughts (Fig. 6), 
especially for the SW+CP+SE droughts that are partic- 
ularly prevalent during the exceptionally dry twelfth and 
thirteenth centuries. Only the SW+CP+NW combina- 
tion shows a small and marginally significant (p < 0.10) 
trend toward decreased occurrence (—0.47 century -1 ) 
over the course of the millennium. For all the pan- 
continental events, the observed percentages of occurrence 
exceed the 99th percentile from our phase-randomized 
resampling. This strongly suggests that pan-continental 
droughts are unlikely to occur as a random consequence 
of independent regional variability, instead indicating 
that they are likely triggered by coherent forcing or 
dynamics acting across regions. 

Average recurrence intervals for single -region (Fig. 7 a) 
and pan-continental (Fig. 7b) droughts increase with the 
transition from the MCA to LIA/MOD period, consis- 
tent with other indicators indicating a shift toward 
wetter conditions in North America. For single-region 
droughts, recurrence intervals range from roughly 2 to 
3 years; for pan-continental droughts, the recurrence 
intervals are all 10 years or more, except for the most 
common pattern, SW+CP+SE, which recurs on aver- 
age every 7 or 9 years during the MCA and LIA/MOD 
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Regional ACFs 

SW CP 



NW SE 



Fig. 5. ACF for each regional PDSI lime series for the MCA (brown solid lines; 1000-1500 CE) and LIA/MOD 
(green solid lines; 1501-2005 CE) periods: (a) SW, (b) CP, (c) NW, and (d) SE. Significance thresholds (dashed lines) 
are based on twice the large-lag standard error of the ACF. Prior to calculating the ACFs, time series were adjusted to 
a mean of zero to remove the influence of shifts in the average baseline. 


periods, respectively. Four-region (SW+CP+NW+SE) 
droughts are especially rare, occurring on average 1 
(20 yr) _1 over the last 500 years. 

c. Forcing of pan-continental droughts 

Spearman rank correlations between the NADA and 
the climate indices (SOI, PDO, and AMO) offer strong 
evidence for the potentially important role of these 
climate modes in forcing regional and pan-continental 
droughts (Fig. 8). Positive values of the SOI have a 
strong association with drought in the SW and southern 
parts of the CP; the influence of the SOI weakens with 
the transition from winter (Fig. 8a) to spring (Fig. 8b). 
Positive values of the PDO in winter and spring cor- 
relate strongly with drought in the NW and SE and 


pluvial conditions in the SW (Figs. 8c, d). The PDO 
correlation pattern has some similarity to the SOI 
correlation, which is not surprising given that the PDO 
is negatively correlated with the SOI (Newman et al. 
2003). It is worth noting, however, that the major cen- 
ters of strong correlation differ between the two cor- 
relation modes. Specifically, the PDO over western 
North America includes a NW-SW dipole that the 
SOI correlation lacks. Correlations with the AMO 
(Figs. 8e,f) are generally negative and widespread, 
with strongest loadings in the SW, CP, and SE (Figs. 
8e,f), indicating the tendency for a positive AMO 
(i.e., warm Atlantic SSTs) to favor widespread drought 
(Kushnir et al. 2010; McCabe et al. 2004; Nigam et al. 
2011 ). 
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Fig. 6. Number of pan-continental drought years (PDSI £ —0.5) in each century, calculated 
for all possible pan-continental drought patterns. To ensure adequate sampling of events, 
drought years were allowed to overlap between the three- and four-region drought categories. 
The total percent occurrences of each drought pattern, calculated over the full time period 
available (1000-2005 CE), are indicated in the figure legend. 
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The ACFs of the climate indices (Fig. 9) offer some 
clues regarding differences in drought variability (Fig. 3) 
and drought persistence (Fig. 5) between the SW and 
CP regions. Autocorrelation in the SOI, which loads 
strongly over the SW (Fig. 8a), varies between signifi- 
cantly negative (lags at 2 and 7 years) and significantly 
positive (5 years) values, which is an autocorrelation 
structure typical of oscillatory phenomena. In other 
words, any SOI-forced drought in the SW is likely to be 
followed by an opposite-sign SOI event within 2 years, 
bringing relatively wet conditions that oppose the initial 
drought, limiting drought persistence. In contrast, the 
AMO has significant persistence that slowly decays out 
to 5 years, with little evidence for any strongly oscilla- 
tory behavior. Given the apparent importance of the 


AMO for drought in the CP (Figs. 8e,f), this suggests 
that this SST pattern may contribute to the exceptional 
persistence of droughts over the CP (persistence in the 
CP is higher than any of the other three regions, espe- 
cially during the MCA). 

For the single-region drought composites, from 30 to 
almost 60 drought-event years are available during the 
overlapping time period between the NADA and the 
climate indices (Table 1). Composite climate index 
values for single-region droughts are shown in Fig. 10; 
composites passing the 90th and 95th percentile esti- 
mated confidence limits are marked with single and 
double asterisks, respectively. Composited SOI values 
are significant and positive for the SW, showing a stronger 
winter signal consistent with the correlation map of Fig. 8. 


Drought Recurrence Intervals 


Single Region Droughts 




Fig. 7. Drought recurrence intervals for (a) single-region and (b) pan-continental droughts, calculated separately 
for the MCA and LIA/MOD intervals. Note that the ranges of the y axes are different in (a) and (b) because of the 
much larger recurrence intervals associated with the pan-continental droughts. 




392 


JOURNAL OF CLIMATE 


Volume 27 


Teleconnections (p) 


SOI (DJF) SOI (MAM) 





-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 

Fig. 8. Correlations between gridpoint PDSI from the NADA and: SOI in (a) DJF and (b) MAM, PDO in (c) DJF 
and (d) MAM, and AMO in (e) MAM and (f) JJA. For the PDO and AMO, correlations are based on filtered (10-yr 
LOWESS) PDSI and climate indices. 


PDO composites are positive and significant for droughts 
in the CP, NW, and the SE during DJF. The MAM PDO 
signal is not significant, and this may be due to the large 
number of drought years in the SE during the nineteenth 
century, reflecting the long-term wetting trend that is 
unlikely to be related to any changes in the PDO. The 
AMO composites show significant positive index values 
for the CP and SE in both the spring and summer, but no 
significant association with SW droughts, despite the 
strong negative correlation in Fig. 8 and previous 


research identifying a strong AMO-drought link in this 
region (McCabe et al. 2004, 2008). 

Fewer composite years are available for pan-continental 
droughts (Table 1), reflecting their relative scarcity, es- 
pecially during the instrumental period. Three drought 
patterns appear to have significant climate index compos- 
ites (Fig. 11): SW+CP+SE (+SOI in DJF), SW+CP+NW 
(+SOI in DJF and +PDO in DJF), and CP+NW+SE 
(+PDO in DJF; +PDO in MAM; +AMO in MAM, 
and +AMO in JJA). Interestingly, these are the most 
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Fig. 9. ACF for: (a) SOI (1884-2012), (b) PDO (1900-2012), and (c) AMO (1854-2012). Significance thresholds (dashed lines) are 
based on 2 times the large -lag standard error of the autocorrelation function. Brown lines are the function based on the MAM (SOI and 
PDO) and JJA (AMO) seasonal-average time series; green lines are for DJF (SOI and PDO) and MAM (AMO). Prior to calculating the 
ACFs, time series were adjusted to a mean of zero to remove the influence of shifts in the average baseline. 


common pan-continental drought patterns to occur in 
the last thousand years. This suggests that there may be 
some predictability for these types of pan-continental 
droughts at leads of one season or more, given that the 
current states of the AMO and PDO can be estimated 
and the SOI can be predicted a few seasons in advance. 
These teleconnection patterns could also explain the 
enhanced occurrence of these droughts during the MCA 
(e.g., the 1100s; Fig. 6), consistent with previous work 
suggesting possible SST forcing of megadrought activity 
(Burgman et al. 2010; Feng et al. 2008; Graham et al. 
2011; Seager et al. 2007a, 2008; Oglesby et al. 2012). 

4. Discussion and conclusions 

Droughts are one of the most expensive natural di- 
sasters to affect the United States, causing significant 
harm to agriculture, human health, and the economy at 
large (Ross and Lott 2003). The 2011 drought, a regional 


event centered in the Southwest and southern plains, 
cost $12 billion and contributed to 95 deaths (NCDC 
2013c). The impact of the 2012 drought, a pan-continental 
event, was much more severe. During 2012, 80% of the 
agricultural land in the United States experienced 
drought conditions, the largest area in a single year since 
the 1950s (USDA 2013), and corn yields were 26% be- 
low expectations (Hoerling et al. 2013). The drought 
also contributed toward making 2012 the third all-time 
highest record year for area burned (NCDC 2013b). 
Ultimately, the 2012 drought was the second most ex- 
pensive natural disaster of 2012 (after Hurricane Sandy), 
causing $30 billion dollars in damages and 123 deaths 
(NCDC 2013c), underlining the significant stresses and 
management challenges imposed by these pan-continental 
events. 

Much effort has been devoted to understanding the 
causes of regional droughts in the past (e.g., Coats et al. 
2013b; Cook et al. 2009; Ellis et al. 2010; Herweijer et al. 


Teleconnection Composites (Single-Region Droughts) 


SOI Composites PDO Composites AMO Composites 



Fig. 10. Composite averages of climate indices during drought years (PDSI £ —0.5) for each region: (a) SOI, (b) PDO, and (c) AMO. 
Significant composites, based on the resampling described in section 2, are indicated by single (90th) and double (95th) asterisks. 
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Fig. 11. As in Fig. 10, but for the pan-continental drought patterns. 


2006; Hoerling et al. 2009; Schubert et al. 2004b) and 
how they may change with increased greenhouse gas 
forcing in the future (Seager et al. 2007b; Seager and 
Vecchi 2010; Seager et al. 2013). Few studies, however, 
have focused on the occurrence and variability of multi- 
regional, pan-continental droughts. We find that pan- 
continental droughts have occurred in various flavors 
over the last millennium, and that the 2012 drought is 
similar to the most common type of pan-continental 
drought pattern found in the paleo record (SW+CP+SE). 
We also find that the occurrences of certain pan-continental 
drought patterns (SW + CP + SE, SW + CP+NW, and 
CP+NW+SE) are strongly linked to well-known tele- 
connections originating from the Pacific and Atlantic 
Oceans. Teleconnections are therefore capable of influ- 
encing the probability of occurrence of pan-continental 
droughts, although any particular event could also be 
influenced by internal atmospheric variability unrelated 
to ocean conditions. For example, Hoerling et al. (2013) 
concluded that the 2012 drought was made more likely 
by slowly varying SSTs, but that the particulars of the 
summer 2012 drought in the Plains was strongly con- 
trolled by seemingly random circulation anomalies. In- 
terestingly, our analyses of the CP and SW also found 
that the MCA megadroughts were generally not syn- 
chronous in time across these two regions and also ex- 
hibited significantly different persistence characteristics. 
This result is not necessarily surprising, given the very 
different climatologies and drought teleconnections as- 
sociated with each region. Our analysis does indicate, 
however, that the megadroughts in North America should 
not be viewed as a single, coherent phenomenon af- 
fecting the entire western United States, but rather as 
a superposition of potentially separable regional events. 

Understanding the underlying dynamics and predict- 
ability of pan-continental droughts remains a significant 
challenge. Interpretations from our teleconnection anal- 
ysis are limited by the relative paucity of pan-continental 


droughts available from the observational era. For ex- 
ample, only six SW+CP+NW+SE droughts overlap 
with the PDO record (Table 1), making it difficult to 
determine whether the full range of climate conditions 
and teleconnections associated with these droughts is 
being fully captured. This limitation could be potentially 
addressed through the use of proxy reconstructions 
of these major teleconnection modes (e.g., Cobb et al. 
2003, 2013; D’Arrigo et al. 2001; Emilc-Geay et al. 
2013a, b; Gray et al. 2004; Li et al. 2011), allowing for 
greater temporal overlap with the NAD A. It is difficult, 
however, to find reconstructions that are continuous in 
time, precisely dated enough to allow for comparison 
with the NADA, and (most critically) based on com- 
pletely independent proxy records. As an alternative, 
general circulation model simulations, such as those 
organized under phase 5 of the Coupled Model Inter- 
comparison Project (CMIP5; Taylor et al. 2012), may 
provide another opportunity to investigate the dynamics 
and teleconnections driving pan-continental droughts. 
Nevertheless, the efficacy of these models to address the 
pan-continental drought problem depends crucially on 
their ability to skillfully reproduce the proper telecon- 
nections, feedbacks, and time scales of response (e.g.. 
Coats et al. 2013a). 

Given the ubiquity of pan-continental droughts over 
the last millennium, it is likely that they will be an im- 
portant and significant feature of North American hydro- 
climate in the decades and centuries to come, with potential 
changes in their severity and/or occurrence with in- 
creased greenhouse gas forcing, southwestern North 
America and the plains, for example, are expected to 
continue drying into the future (Seager et al. 2007b, 
2013), with greenhouse warming causing both precipita- 
tion declines and increased evaporative demand in these 
regions. From a pan-continental perspective, warmer 
temperatures may be especially important, as tempera- 
ture responses to increased greenhouse gases are likely 
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to be more widespread and homogenous than pre- 
cipitation changes (e.g., Knutti and Sedlacek 2013). This 
spatially extensive warming, and the accompanying in- 
crease in evaporative demand, represents a broad forc- 
ing that may increase the spatial coherence of droughts 
across regions, increasing the likelihood that multiple 
regions will experience drought simultaneously and 
amplifying the impacts on ecosystems, agriculture, and 
society at large, as was seen in 2012. Nevertheless, a 
complete quantification of the impacts of greenhouse 
warming on pan-continental droughts will require the 
analysis of climate model simulations that can accu- 
rately reproduce the underlying dynamics governing 
drought variability over North America. 
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